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Abstract. We have studied the critical behaviour of a doped Mott insulator near the metal-insulator 
transition for the infinite-dimensional Hubbard model using a linearized form of dynamical mean-field 
theory. The discontinuity in the chemical potential in the change from hole to electron doping, for U 
larger than a critical value Uc, has been calculated analytically and is found to be in good agreement with 
the results of numerical methods. We have also derived analytic expressions for the compressibility, the 
quasiparticle weight, the double occupancy and the local spin susceptibility near half-filling as functions 
of the on-site Coulomb interaction and the doping. 

PACS. 71.10.Fd Lattice fermion models (Hubbard model, etc.) - 71.27.-|-a Strongly correlated electron 
systems; heavy fermions - 71.30.-|-h Metal-insulator transitions and other electronic transitions 

1 Introduction namical mean-field theory (DMFT) |Q to generic micro- 

scopic models, such as the single-band Hubbard model. 

mi TT 1 1 1 i 1 • 1 i i -i- /n TTrn\ • lu thls aporoach the lattice problem is mapped onto an 

The Mott-Hubbard metal-msulator transition (MIT) is ^ 

i 1 i i-£i2 li 11 11 1-11 impurity problem where a correlated impurity site is em- 

an important but dirncuit many-body problem wnicn has 

, ii-ii •iOiri--j2 i 1 1 bedded in an effective uncorrelated medium that has to be 

been studied extensively |ll| . bignihcant progress has been 

, . J . , , 1- ,• r J determined self-consistently. The mapping can be shown 

achieved m recent years through the application of dy- 
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to be exact for certain models with local interactions in are in very good agreement with most accurate numerical 

the limit of infinite dimensions and constitutes an ap- estimates. 

proximate method for finite-dimensional models. Several Here we show that the linearized DMFT can be ex- 
techniques have been applied to deal with the effective im- tended to the non-symmetric case. This allows for a com- 
purity problem, including iterated perturbation theory |^ , prehensive analytical investigation of the critical behaviour 
non-crossing approximation projective self-consistent as a function of doping. In particular, we have calculated 
method (PSCM) 1^, quantum Monte Carlo (QMC) Q, the discontinuity in the chemical potential on changing 
exact diagonalization (ED) 0| and numerical renormaliza- from hole to electron doping. The analytic result is in 
tion group (NRG) |^,^ techniques. These methods have good agreement with the results of numerical methods, 
different strengths and limitations, and apparent discrep- NRG as well as ED. It also agrees well with the result of 
ancies between some of the results have led to a contro- PSCM ||]. In addition, we have analytically determined 
versial discussion on the nature of the MIT in the infinite- the compressibility, the quasiparticle weight, the double 
dimensional Hubbard model at low temperatures ||l0|,^, occupancy and the local spin susceptibility near half- filling 



|12|,|13, 14|. Recently, however, most of these differences have as functions of the on-site Coulomb interaction and the 



been resolved [|5|. doping. These are difficult to calculate using the numeri- 

TT .11 J- fj-T, 4- v u u • cal methods mentioned above. 

However, not all aspects of the transition have been in- 
vestigated thoroughly. In particular, the critical behaviour 

of the doped Mott insulator has not received so much at- 2 Linearized dynamical mean-field theory 

tention as it is more difficult to solve the model off half- 
filling using purely numerical methods. Iterated perturba- 
tion theory, which provides a reasonable description of the 
MIT at half- filling |Q, is less suitable for calculations in 
the doped case since ad hoc modifications of the approach 



We consider the single-band Hubbard model at zero tem- 
perature on the Bethe lattice [l9| ] with infinite connectiv- 
ity q i—f oo 



<i,j>.a i 

become necessary away from half-filling ||l6). The nearest-neighbor hopping is scaled as usual: Uj = 

Recently, we have shown 0[l| that a linearized form §' the single-band Hubbard model on the Bethe lat- 
of the DMFT provides a simple and attractive technique ^ice, the DMFT self-consistency equation is simply given 



to obtain approximate but analytical results for the crit- 
ical regime. Using this approach for the single- and for 
the two-band Hubbard model, it has been demonstrated where G{z) is the local Green's function of the Hubbard 



by i 

go{z)-^ = z + fi-t^G{z), (2) 



that the predictions for the critical coupling of the MIT model and (Jo (2^) (^~£/~^(-2^)) ^ is the non-interacting 
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Green's function of an effective single-impurity Ander- bandwidth Anderson model [ pT[ 
son model defined by the hybridization function /^(z). 

The DMFT self-consistency cycle starts with a guess for o' 

_ +e,y2cic„ + vy2ific^+ciu) (s) 

the hybridization function |E0|. The solution of the cor- ^ ^ 

responding impurity model yields the impurity Green's with — and £/ = — ^. The hybridization strength 

function G{z) to be identified with the on-site Green's V has be determined from the self-consistency equation 

function of the lattice model. From the self-consistency which takes the simple form 
equation (0) Go{z) and thus a new hybridization function 

fZ = V^. (4) 

A{z) can be calculated. The cycles have to be repeated 

until self-consistency is achieved. For a more detailed discussion and for an extension of the 

approach to the two-band model see refs. [ p7| , p^ ]. 

The application of numerical methods for solving the Nqw we calculate the critical value of U for the MIT at 

effective impurity problem (see ref. Ill, for example) have half-filling. In this case, the chemical potential is fixed to 

established that at moderate coupling the density of states ^ ^ !l due to particle-hole symmetry. Close to the MIT 

in the metallic phase is characterized by a three-peak y ig gmall, and F ^ for U ^ Uc- The weight Z can 

structure consisting of the lower and the upper Hubbard be calculated from the impurity spectral function and is 

bands and a quasiparticle resonance near the Fermi en- given by j2l| (see also Appendix ^ 
ergy in addition. The quasiparticle peak is more or less 

(5) 

isolated from the Hubbard bands. On approaching the 

Mott-insulating state with increasing interaction strength, up to second order in V. From eqs. (|,|), we obtain the 

the quasiparticle peak becomes extremely narrow and fi- critical value 
nally vanishes for U = Ur- 

Uc = 6t. (6) 

The linearized dynamical mean- field theory fo- This result is in good agreement with the best numerical 

cuses on this critical regime close to the transition. Here, estimates of NRG {Uc = 5.88^) jl3| and ED (Uc — 5.87t) 

the quasiparticle peak is approximated by a single pole as well as PSGM {Uc — 5.84i) [Q. 

at the Fermi level, i. e. G{z) — ^ near the Fermi level When we solve the self-consistency equation ^ by it- 

with a small weight Z ^ as [/ — > [/c. Correspond- eration, T^^ increases exponentially with iteration number 

ingly, the hybridization function A{z) is a one-pole func- for U < Uc. In this case, the single-pole approximation 

tion A{z) — This effectively represents an approxi- for A{z) breaks down. On the other hand, for U > Uc, 

mate mapping of the model (|l|) onto a two-site or zero- the hybridization strength decreases exponentially to 
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give the self-consistent value — corresponding to the up to second order in V. We also calculate the static local 



insulating solution. 

3 Critical behaviour of the single-band 
Hubbard model 

3.1 Physical quantities 

For more extensive studies of the critical region of the 



spin susceptibility, xo = Xm(0)j with 



(13) 



where Si are the spin raising and lowering operators at 
the same site i and where we set (17/15)^/2 = 1. At T = 0, 
Xo is given by 

m,\S+\Eo)\' + \{E^\S-\Eo)\' 



En — Eq 



(14) 



Mott MIT, we need the result for the quasiparticle weight where \En) are the eigenstates with the eigenenergies En, 

and \Eo) is the ground state. By using the results for the 
2-site Anderson model in eq.(|lj), we obtain (see Appendix 



up to fourth order in V (see Appendix 

Z = V'F{U,fi)-V^G{U,fi), 
where 



(7) 



Xo 



2^2 



F{u,^l) = — 



oq 



2^2 n{U-fi) 2([/-^)2' 
24 22 



2^4 ^3([/_^) ^2([/_^)2 

24 29 



(8) 



(9) 



(15) 



We now invoke the self-consistency equation (Q) and 
substitute from eq.@ into eqs.(0,|l^,|l|,|l5|). This yields 



To determine the doping dependence, we also calculate 
the impurity occupation number up to second order in V 
(see Appendix ^) 



z{u,^i) = [F{u,^i) 



n = l + V^D{U,ti), 



(10) 



where 



<2 ; t^G{U, fi) ' 
n{U,fi) = 1 + Z(U,fi)fD{U,fi), 

d{U,ti)^2eZ{U.ji)/{U-ii)\ 

Xo{U,y) = ti{U-y)/{2Ut^Z{U,tj)). 

3.2 Critical behaviour at half-filling 



(16) 
(17) 
(18) 
(19) 



^(f/,/i) 



(11) 



First we discuss the critical behaviour near the Mott MIT 
at half- filling [f?!. The chemical potential is fixed to /i = 



We note that, for ^ | -j, D | and thus n | 1. The ^ due to particle- hole symmetry. In fact, cq.(n^) yields 



double occupancy d = {fi^fii) is given by (see Appendix n(J7, y) = 1 for aU U as £>([/, -j) = from eq.(|Tl|). Sub 



0) 



stituting = Y eqs. (|16| , |18| , |19D , we obtain 



2V' 



(12) 



18 



11 



u 



Uc 



(20) 
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f^^jj-j — ^ f I j (21) ^ discontinuity at ?i = 1: At half-filling the system is a 

11 \ Uc ) 

, , / TT \~'^ Mott insulator for /i_ < /i < /i+ where /i_ = /i_ (?7), /i+ = 

Xo{U) = (1- — ) , (22) 

24< V UcJ /^+(C/), and ^^{U) < fi+{U) for U > Uc. The jump A^i = 

near [/, for U < C/,. The result for Z, eq.®, has already _ ^_ ^an be calculated within the linearized DMFT 

been obtained by Bulla and Potthoff 0. Using eq.® exploiting the fact that the MIT is characterized by a 

the compressibility, k = (f^)„, is given by vanishing quasiparticle weight. Setting Z = in eq.®. 



.(f/)=t^Z(C/)^^(^'^) 



(23) 



we obtain an equation to determine the MIT point: 



at half-filling for U < Uc- Substituting eqs.(|l|,|2^) into FiU,n) = ^. (25) 



eq-(|23|), we obtain ^ p^^^ plotted as a function of for several 

^([/) = M _ ^1 ^ (24) values of t/. When U < Uc, F{U, fi) > ^ for aU ^, and the 

33i Y Uc J 

system is metallic. When U > Uc, F{U, jj) — ^ for IJ. ^ fJ.± 



near Uc for U < Uc- 

Due to local charge fluctuations the double occupancy 
is finite even in the insulating phase with d{U) —* only 
for U oo. The linearized DMFT, however, predicts a 
vanishing double occupancy for U Uc. This indicates 
that the approach misses these local fluctuations. We nev- 
ertheless believe that the trends in the metallic phase are 
correct on the mean-field level. This situation may be anal- 
ogous to the mean-field theory of ferromagnetism in the 
Heisenberg model, where the fluctuations in the magneti- 
zation and short-range order cannot be taken into account 



> /i-), and the system is insulating for /i_ < /i < 
while it is metallic for /i < and /i > Wc notice that, 
for U — Uc, F{U, ii) as a function of fi is at its minimum 
for A* = 3" consistent with the considerations in ref. fl^ . 

Using eq.(§) we can solve eq.(^) to get the explicit U 
dependence of fi± for U > Uc 




where u = jj- > 1. For the discontinuous jump Aji = 



We also note that the local spin susceptibility diverges in a,, rr i i ^ ^ . / i f ^ \ 1 fo'7\ 



the limit of the MIT U Uc. Even in this limit, however, 
the uniform spin susceptibility is finite due to the effect of 
the super-exchange interaction J = 2i^/C/ [p[. 




Note that we have n ^ 1 with n < 1 (n > 1) from eq.(jl^) 
for fi fi^ {fi — !■ /i_|_) from the metallic side. For U > Uc 
close to Uc, eq.(|27|) yields 



3.3 Discontinuity in the chemical potential An = -^—Ur\ — — 1 (28) 

V38 M Uc ' 

For U < Uc the chemical potential fi as a function of n is In Fig. |^ we have plotted Afi, as given by eq.(|2^), as a 
continuous at rt = 1. For U > Uc, on the other hand, /i(n) function of U. 
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We have also calculated Afj, numerically by using the within the quasiparticle peak which in turn lies within 
NRG and the ED method. We solved the full DMFT equa- the Mott-Hubbard gap. The same applies to the slightly 
tion numerically and obtained the occupation number n as electron-doped case. This interpretation is confirmed by 
a function of /i. When /i approaches /i+ (/i-) from above more extensive NRG and ED calculations, 
(below), n approaches unity from above (below) and n — 1 

for < ^ < /!+. We also confirmed that, at /i = /i±, the 3.4 Doping dependence 

groundstate changes from singlet (/i < /i_ or /i > /i+) to 

Finally, we discuss the critical behaviour of the doped 

doublet (/i- < 1-1 < The ED calculations were done 

metallic system for U > Uc close to the Mott insulat- 

for finite cluster sizes rig up to = 11. An extrapolation 

ing phase. Substituting eq.(p6[) into eq.(|l7|), we obtain the 

of the results yields the ti^ — > oo extrapolated value of Afi 

doping dependence of the quasiparticle weight for n 1 

as shown in the inset of Fig. p|. 

and [/ > Uc 

As can be seen from Fig. ||, the result from the lin- 

Z{U,n) ^ C{U)\n-l\, (29) 

earized DMFT is in good agreement with the results of 

the NRG and the ED. It also agrees well with the result where the coefHcient C(U) is given by 

from the PSCM H. We note that the result for ^± from r^nT\-i n/rr m 00^2 /on\ 

the linearized DMFT seems to be less satisfying for the 

Due to particle-hole symmetry we obtain the same result 

strong-coupling limit U ^ Uc- This may be explained by 

. , 1 1 r 1 1 for both, n < 1 and n > 1. 

the fact that here fi± is very close to the edge of the lower 

, , , , , . , , , , rr r Eliminating Z from eqs. (flq,0) , the occupation mim- 

or upper Hubbard band, respectively, and that the etrect 01 ^ — " — ' 

, , , ., 1 1.1- , , ■ , ,. . iT^^TT^m ber can be expressed as a function of [/ and u. For [/ > Uc 
the bandwidth, which is neglected m the linearized DMl" i 

J— 3, , . J— I and close to the MIT we have 

| |17[ |, becomes important 

In the case of half-filling, as U approaches Uc from be- 

n{U, m) = 1 + (^F{U, /i) - c!(u]^J) - ^^^^ 

low, we have the picture of well separated lower and upper r^,. , • u j-t, -un-i- • j-t, r v 1 

' ^ ihis also yields the compressibility m the limit n ^ 1 tor 

Hubbard bands with a quasiparticle peak in between. The 

U > Uc 

quasiparticle peak disappears at U — Uc leaving an insu- 



k{U} - 



(32) 



lator with a finite gap. For U > Uc, Fisher et al. Q give G{U,^±) dfi 

a similar picture for the slightly doped case in terms of Again, due to particle- hole symmetry, the result eq.(|3 

quasiparticle in-gap states. The jump zi/i from hole to for k is the same for both, n < 1 and n > 1. 
electron doping tends to zero as U Uc from above, so Substituting eqs.( p6| , p9| ) into eq.(p^, the double occu- 

the Fermi level in the slightly hole-doped case must lie pancy is obtained as a function of U and n for n I and 
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Fig. 1. The ^ dependence of F{U, /i) defined in eq.M for several values of U. The critical interaction is Uc = 6t. We set t — 1. 



Aju 




Fig. 2. The discontinuity of the chemical potential = ^+ — as a function of U for U > Uc, as obtained from the 
linearized DMFT (solid line), from the exact diagonalization method (closed circles), from the projective self-consistent method 
(open squares) jsj and from the numerical renormalization group method (crosses). Inset shows an extrapolation of the ED 
results with system sizes Us — 6, 8, 10 and 11. We set t — 1. 



Y. Ono, R. Bulla. A. C. Hewson, M. PotthoflF: Critical behaviour near the Mott metal-insulator transition 



d{U,n)=a±iU)\n-l\, (33) p^u)^^J—-l. (40) 



approach |23,24]. However, the coefficients are different. 



4 Conclusion and discussion 



where The critical properties near Uc as resulting from the lin- 

16{U — ii±)^UAfx earized DMFT are similar to those of the Brinkman-Rice 

and where + and — stand for n > 1 and n < 1, respec- 
tively. For a given interaction strength U we find a-|_ > a_ . 
This means that when doping the system away from half- 
filling, the increase of the double occupancy for n > 1 is 
stronger than its decrease for n < 1. As mentioned be- 
fore, there is a finite double occupancy even for n = 1 
and U > Uc in contrast to the result eq.(|3^) from the 
linearized DMFT. Nevertheless, we believe that the trend 
in the metallic state, i.e. the coefficient a±, is physically 
significant. 

Substituting eqs.( |2^j29| ) into eg. (p^ , the local spin sus- 
ceptibility is obtained as a function of U and n for n — )■ 1 
and U > Uc 



We have shown that analytical calculations of the critical 
behaviour in the parameter region close to the Mott MIT 
can be carried out by using a simple linearized version of 
the DMFT. This approach gives rather reliable estimates 
for the critical parameters at the transition point. Consid- 
ering the single-band Hubbard model, we have calculated 
the discontinuity in the chemical potential /i for the change 
from hole to electron doping. The analytical result is in 
good agreement with the numerical estimates from the 
PSCM, the NRG and the ED. The resuhs for smaU hole 
and electron doping and for U > Uc can be interpreted in 
terms of quasiparticle in-gap states which lie within the 

1 Mott- Hubbard gap. We have also made predictions for the 

where ° ^ 



Xo(U,n)=(3{U)\n-l\-\ (35) 



compressibility, the quasiparticle weight, the double occu- 

/?([/) ^ "^"^^ (36) 

C/^ — (/i/i) pancy and the local spin susceptibility near half-filling as 

In Fig. I C-\ K, a^^ and /?, as given by eqs.(|^,||, functions of U and the doping. It would be interesting to 

||,|^), are plotted as functions of U. To see the critical see whether these predictions can be confirmed by some of 

properties near U^ we substitute eq. (||) into eqs. (|o|j3||3|, t^^e more comprehensive numerical reahzations of DMFT, 

^), and obtain for U > Uc whether there is any experimental evidence in support 

Ifi fu of them (2|. 

C(C/)-i = (37) 

3\/38 V f^c The simplicity of our approach should also enable sim- 

<U)^^(^-l], (38) 



33t \Uc / ' ^^^^ calculations to be made for more general models, such 
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Fig. 3. The compressibility k, the coefRcient C of the quasi 
occupancy d = a±\n — 1| with ± for n < 1 and the coefficient f: 
as functions of U for U > Uc ~ 6t obtained from the linearized 

as multi-band Hubbard- type models, where we can have 
both Mott-Hubbard type and charge-transfer type metal- 
insulator transitions [ p6| . Even for more complex multi- 
band models one might obtain analytical expressions as 
has been demonstrated recently for the five critical pa- 
rameters in the two-band Hubbard model . This makes 
the linearized DMFT a powerful technique for explor- 
ing general trends and phase diagrams within a high- 
dimensional parameter space. An interesting question to 
examine would be the role of in-gap states in the two-band 
Hubbard model (equivalent to the d-p model) near the 
MIT. In this parameter regime the results of the model 
should be relevant for a description of the behaviour of 
weakly doped high-Tc materials. 



aarticle weight Z = C\n — 1\, the coefficient a± of the double 
' of the local spin susceptibility xo = l3\n — near half-filling 
DMFT. We set t = 1. 
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Appendix A 

Here we discuss the two-site Anderson model eq.(||) in 
the Umit V ~^ 0. We assume that the conduction level is 
between the atomic /-level and the upper Hubbard level 
Ef < Ec < Sf + U. Then we define 



A = ec-ef> 0, 

W = ef + U -ec>0. 

The one-electron eigenstates 

|i?±)=a±/+|0)+5±c+|0), 

correspond to the eigenenergies 



(41) 
(42) 



(43) 



£c + e/±x/(£c-£/)2+4y^ 



For small hybridization strength ^ A, eq.(^ 
plified as 



IS Sim- 



E. 



E. 



to fourth order in V, with the corresponding eigenstates 



\E^)=alf+ 



1 - 



A^ 



V 



1 



|0), 



Z\2 



(47) 
(48) 



where 



2 



(49) 



Similarly, we obtain the three-electron (one-hole) eigenen- 
ergies 



E± = ]^[iec + 2,ef + U± ^/{ef + U- Sc)^ + 4V^ )(50) 



For small hybridization strength V'^ <C W, eq.(|50|) is sim- 
plified as 



E+ =ec + 2ef + U 
E_ =2ec + ef 



W W3' 



(51) 
(52) 



to fourth order in V. The corresponding eigenstates are 



V 



W 



^+) ( l--^)/- + c4l4>, (53) 



V 



E^)^aiU + —[l- — ]cA\A) 



W 



M/2 



where 



2 F2 y4 



(54) 



(55) 



and |4) = /+/+C+C+IO). 

The two electron states can be classified as singlets 

(44) or triplets. In the triplet state, the spatial part of the 
wavefunction is antisymmetric and the interaction U has 
no effect. Then the total energy of the triplet state is given 
by £'+ + E^ = Ec + Ef. There are three possible singlet 

(45) states which can be written by the linear combination of 

(46) the states, 



i</>i)-^(4./r-4/T')io>> 

\h)^c+c+\0), 

\h) = f^ftm- 



(56) 
(57) 
(58) 



The eigenenergies are given by the solutions of the equa- 
tion, 

E-Ec-Ef -V2V -V2V 
-V2V E~2ec 
-V2V E-2Ef-U 

To fourth order in V , the ground state eigenenergy is 



0. 



(59) 



1 1 



Eo=E. + Ej-2V^[- + -][l- — - — ]m 



2V^ 2V' 



Z\2 



M/2 
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and the corresponding singlet ground state is 



X 



\Eo) = ao { - 



1- 



W 



1- 



aw) 



17_ 



24 



22 



+ 



24 



A^W A^W^ AW^ 



AW 



(61) 



\{E-\MEo)f 



wy 

(65) 



with 



Z\4 



A^ 
+ 



AW 



16 



+ 



4 

' W^ 
22 



+ ■ 



Yl 

2 

32 



+ 



24 



— —\ 



= W2, 



A^W A^W^ AW^ ' 

(66) 



3 
W^ 



(62) 



^A^ A^W A'^W^ AW^ 

Now we calculate the /-electron Green's function of 

this model. When a f , t electron is removed from the , -ii • i x- i -j. ^- -^i 

•' there are two possible smgle-particle excitations with 



to fourth order in V. 

When a /, t electron is added to the ground state \E\), 
possible final states are \E-) and |£'+). Correspondingly, 



ex- 



ground state \Eq), there are two possible final states: \Ej^) 
and \E-). Correspondingly, there are two possible single- 
hole excitations with excitation energies, 



5 



3 2 
A^W 



+ 



A^ A'^W AW^ 



4 



(63) 



E_-Eo = -£c + 



1 2 
A'^W 



= -£2, 



3 

^3 



A^W AW^ 



4 



(64) 



to fourth order in V. The matrix elements for these tran- 
sitions are 

2V^ 



{E+\f^\Eo} = 



{E_\f^\Eo) = 



aao 



71" 

w 



AW 

Yl 

V 
A 

2y2 
^ w^ 



1 



_ ^ 
" w^ 

Yl 

2T/2 \ 
Aw) 



2^2 
AW 



which yield the transition probabilities: 



\{E+\MEo) 



4 2 



Z\2 AW W^ 



2 



citation energies, 



E. 



En 



2 1 
4 



Z\3 /\2^ /\^2 



3 



£3, 



(67) 



E+- Eo = £f + U + V'^ 

= £4, 



2 3 



/\3 ^^2 



5 



(68) 



to fourth order in V. The matrix elements for these tran- 
sitions are 



{E-\h\Eo) 



y/2\w\ W^ 

^ ^ V ~^^Aw) y 



{EMEo) = ^{l 



X 1 



_ ^ 
AW 

Yl 
w^ 



2T/2 



which yield the transition probabilities: 



\{E-\MEo)f 



T/2 / 4 



2 ' ' W"2") 



+ 



+ 



32 



22 



2^2 \ 
AW J 



Yl 

2 
16 



24 

^4 zi2^2 + AW^ VF-i y 
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= W3, (69) 2. For a review, see A. Georges, G. Kothar, W. Krauth, M. J. 

I(^+I/tI^o>P = 



1 f 2 4 3 \ Rozenberg, Rev. Mod. Phys. 68, 13 (1996). 



2 2 yA^ AW J 2 

^2 24 22 24 17 \ ^- ^etzner, D. VoUhardt, Phys. Rev. Lett. 62, 324 (1989). 



Z\4 A^W A^W^ AW^ W^y 4. Th. Pruschke, M. Jarrell, J.K. Freericks, Adv. Phys. 44, 187 

= W4, (70) 

^ ' (1995). 

5. D. S. Fisher, G. Kotliar, G. Moeller, Phys. Rev. B52, 17112 
to lourtn order m V . i i j j : 



(1995). 

6. M. Jarrell, Phys. Rev. Lett. 69, 168 (1992). 

7. M. Caffarel, W. Krauth, Phys. Rev. Lett. 72, 1545 (1994). 
Ga{z) = X! ^ 1\. ■ C^-^) 8. O. Sakai, Y. Kuramoto, Solid State Commun. 89, 307 

(1994). 

9. R. Bulla, A.C. Hewson, Th. Pruschke, J. Phys. Cond. Mat. 
10, 8365 (1998). 

10. J. Schlipf, M. Jarrell, P.G.J, van Dongen, N. Bliimer, S. 
Kehrein, Th. Pruschke, D. VoUhardt, Phys. Rev. Lett. 82, 



From eqs.( p3[ - [70| ), we obtain the /-electron Green's 
function which has four poles: 

4 

Wi 

. z- 

In the limit V 0, high-energy poles at ei « £/ and 
£4 « £/ -f t/ have large residues wi ~ W2 ~ ^, while low- 
energy poles merge together at £2 w £3 ~ with small 
total weight Z = 1112 + w^: 



Z = V'^(— + + 4890 (1999). 

\2A^ AW 2W^) 



/ 29 24 22 24 29 \ J- Rozenberg, R. Chitra, G. Kotliar, Phys. Rev. Lett. 

V 2^4 ^ A^W ^ A^W^ ^ AW^ ^ 2W^ I ' 

^ ^ 83,3498 (1999). 



to fourth order in V. The number of / electrons, n = 12. R. M. Noack, F. Gebhard, Phys. Rev. Lett. 82, 1915 
2{wi + W2), is given by (1999). 



2V'^ ( — 



^ s 13. R. Bulla, Phys. Rev. Lett. 83, 136 (1999). 

) 14. W. Krauth, Phys. Rev. B 62, 6860 (2000). 



3 2 2 3 

( ^ + - 'AW^ ~ W^ ) ' VoUhardt, |sond-mat/0012329 



to fourth order in V . By using eg. (pi]), we obtain the dou- 
ble occupancy, d — {EQ\h f^fi f i\Eq) , to fourth order in V 

2V^ ( / 2 4 6 



16. H. Kajueter, G. Kothar, Phys. Rev. Lett. 77, 131 (1996); 
M. Potthoff, T. Wegner, W. Nolting, Phys. Rev. B 55, 16132 
(1997). 



d = ll-V"^ { — ^ — + — W . 17. R. Bulla, M. Potthoff, Eur. Phys. J. B 13, 257 (2000). 

W^ \ \A^ AW W^ ' ' . V ; 
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